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ABSTRACT 

Context. The Cassini mission in the Satumian system is an outstanding opportunity to improve our knowledge of the satellites of 

Saturn. The data obtained thanks to this mission must be confronted to theoretical models. 

Aims. This paper aims at modeling the rotation of Mimas, with respect to its possible internal structure. 

Methods. For that, we first build different interior models, in considering Mimas as composed of 2 rigid layers with different poros- 
ity. Then we simulate the rotational behavior of these models in a 3-degree of freedom numerical code, in considering complete 
ephemerides of a Mimas whose rotation is disturbed by Saturn. We also estimate the deviation of its longitudinal orientation due to 
tides. 

Results. We expect a signature of the internal structure up to 0.53° in the longitudinal librations and an obliquity between 2 and 3 
arcmin, the exact values depending on the interior. 

Conclusions. The longitudinal librations should be detectable, inverting them to get clues on the internal structure of Mimas is 
challenging. 
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' 1. Introduction 

[ The Cassini spacecraft gives us the unique opportunity to have 
. accurate set of geodetic data for icy satellites of Saturn as 
' for e xample, the sh ape, the gravitational field, the rotational 
, data dlh omasLl20Tol) . The flybys of Mimas have provided high 

■ resolution images of the surface in the finest detail yet seen 
' (iRoatsch et al.L l2009l) . Cassini space c raft h as detected temper- 
. ature inhomogeneities (iHowett et al.L 1201 ih . usually attributed 
' to exogenic process. The theoretical model of Mimas rotational 
] state can be used to interpret the Cassini data and to better un- 
. derstand its interior and evolution. 

' Like for our Moon, Mimas is in synchronous rotation and 
, shows almost the same face towards Saturn. Moreover, it is con- 

■ sidere d to have a large librational amplitude (IComstock & BillsL 
I l2003h . The rotational state of a synchronous body depends on the 
i distribution of mass of the body, and therefore it is a signature 
' of its internal structure. Here, we propose to model the rotation 
, of Mimas considering it as a rigid body. A highly rigid interior 

of Mi mas for most of its history is consiste nt with its un-relaxed 
shape dThomas et al.L 120071: lThomasLl201 Oh . 

Since the distant spacecraft flybys of Mimas do not allow 
the determination of the GM nor the gravity harmonics, the 
mass of Mimas is determined from an analysis of its orbita l 
resonances with Tethys and Methone jJacobson et al.L l2006h . 
Moreover, its internal structure remains uncertain. The mean 
density of 1.15 g cnr^ suggests that Mimas is made of ho- 
mogenous mixture of ice and rocks. The observed shape of 
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Mimas by Voyager has been interpreted as an indication of 
interior mass conc entration which can be eit her due to inter- 
nal differentiation ("Dermott & Thomasl Il988l) or radially vari- 
able porosity ( lEluszkiewica il990l) . However. Cassini observa- 
tions showed that Mimas' shape, although a triaxial ellipsoid, 
is departed slightly from hydrostatic shape and therefore in- 
terpreting the i nterior configu r ation from the shape is limited 
dThomas et al.L I2007L ITh omasL l2010h . In the present study, we 
consider Mimas to be composed of two rigid layers. We consider 
both hydrostatic and non-hydrostatic interior models. The inte- 
rior models co nsidering compa ction of ice-silica particle mix- 
tures dYasui & ArakawaL |2009|) are expe cted to yield realisti c 
principal moment of inertia A < B < C dEluszkiewiczi Il990l) . 
Since Mimas orbits close to its parent planet, the present-day di- 
urnal tidal stresses can be important, and we took the tidal effects 
into account as well. 

The paper is structured as follows: we first model the internal 
structure of Mimas, in considering two different assumptions: 
that Mimas is in hydrostatic equilibrium, and that its ellipsoid 
of gravity is proportional to its ellipsoid of shape. From these 
two assumptions we derive 23 models of Mimas. Then we per- 
form numerical integrations of the rotation of these "Mimasses" 
in a full 3-degree of freedom conservative models. Finally, we 
check the influence of the tides on the equilibrium orientation of 
Mimas' long axis. 

2. Internal structure 

Interior structure models of planets and natural satellites are in 
general non-unique due to the presence of fewer constraints than 
unknowns. For Mimas we have only two constraints: the mean 
radius R and the mean density p, or the Mass m (TabU]). In this 
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Table 1. Physical and dynami cal properties of Mimas, used in the ca lculations. The mea n density p has been calculated from the 
radius of iThomas et al.l (l2007h and the mass of Ijacobson et al.l (l2006l) . A recent paper by iThomasI (|20 1 Ol) slightly shifts the mean 
radius to 198.3, the change has negligible effects. However, we use this last reference for the triaxial shape of Mimas, because the 
rotation, especially the longitudinal motion, is sensitive to the differences between these axes. 



Parameter 


Value 


Source 


Mean motion n 


2435.14429644 rad/y 


TASS1.6 (Vienne & Duriez. 1995) 


Mean radius R 


198.2 km 


(Thomas et al.. 2007) 


Density p 


1150.03 A:g.m-3 


(Thomas et al.. 2007) 


Mass m 


3.7495 X 10" kg 


f Jacobson et al.. 2006) 


Saturn-facing radius a 


207.8 km 


(Thomas, 2010) 


Orbit facing radius b 


196.7 km 


(Thomas, 2010) 


Polar radius c 


190.6 km 


(Thomas, 20 Ij)) 



study where we deal with the rotation, the moment of inertia 
differences are the main point of interest. 

Because we have only few constraints, we prefer to have 
as simple interior structure models as possible. We assume a 
two-layer interior structure model with a rocky core and icy 
mantle. The models with variable porosity suppor ted by recent 
compaction experiment s (lYasui & Arakawa , l2009h . yield realis- 
tic moments of inertia ( lEluszkiewiczl 1 1 99C ). The mean density 
of Mimas is close to the density of water ice and the interior is 
consistent with an icy mantle and small rocky core, alternatively 
Mimas can be homogenous with a variable or constant porosity. 

With Mirnas' figure departed from hydrostatic shape 
(iThomas et al l l2007H Th oma and without the knowledge 

of C22 and J2, we do not know the moment of inertia differ- 
ences. As it is the case for the Moon (see e.g. iLam beck & Pullanl 
(I1980I) ). the internal mass distribution could be a a fossil shape 
which dates back from an earlier orbital position where the tidal 
heating was important with higher orbital eccentricity and obliq- 
uity or a frozen shape following a large impact. Mimas shows a 
heavily cratered surface without signs of geological activity for 
billions of years. Its large free eccentricity gives another reason 
for low internal activity. This anomalously large eccentricity of 
Mimas can be explained by pass age through several past reso- 
nances (iMever & Wisdorril l2008h . Determination of the gravity 
coefficients C22 and J2 - -C20 are necessary to conclude on 
the hydrostatic equilibrium because the figure of Mimas may 
not represent the real flattening or the internal mass determi- 
nat ion as it is the case f or Titan (jZebker, 2009). As suggested 
by Ijohnson et al.l (l2006l) . large impact craters and heating of 
Mimas' hemispheres by Saturn at different amounts may be po- 
tential sources of large non-hydrostatic anomalies that could im- 
pede accurate interpretation of the shape data 



2.1. Hydrostatic approximation 

Since we do not know the gravity coefficients C22 and J2, we will 
use a simple approach. For a two layer interior model the core 
radius Re can be determined if the densities of the rocky core 
and the icy mantle p, are known, i.e. 



Rr^R 



P-Ps 



1/3 



\Pc - Ps 

The moment of Inertia factor (MOI = Ip/{MR^)) is given as: 



(1) 



MOI = - 



(P-P.)"' 



Pi 
P 



(2) 



^\p(Pc-Psf'^ 

The FiglT]shows that in the absence of additional constraints, 
plausible density values of the core and icy shell yield 0.3 < 
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Fig. 1. Variations of the Moment of Inertia factor (MOI) with icy 
shell and rocky core densities. 



MOI < 0.4. The range of MOI is used to estimate the moment 
of inertia of tri-axial Mimas, as shown below. 

For a satellite in hydrostatic equilibrium, MOI is re- 
lated to the fluid Love number kj which describes the re- 
action of the satellite to a p erturbing potential a ft er all 
viscous stresses have rela xed (iMunk & MacDonaldl 1 19601 : 
iHubbard & Anderson! [T978h : 



MOI 



1 +k 



(3) 



and the gravity coefficients C22 and J2 are determined from 
(Rappaport et al..,1997i) : 



C22 - 



Jo 



5kf / ,\ 

—q,- + 0[q-). 



(4) 
(5) 



where - QrR^ l(Gm), Q. being the spin velocity of Mimas, 
equal to its mean motion n since Mimas is in synchronous 
rotation. With the numerical values given in TabHj we have 
qr = 0.01854. 

The differences between the three principal moments of in- 
ertia A < B < C are determined from the definitions of C22 and 
J2, i.e.: 



B-A = AC22MR^, 
C-A ^ (J2 + 2C22)MR\ 
C-B = (72 - 2C22)MR^. 
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The relationship between the mean moment of inertia / 
and the polar moment of inertia C is: 



C 



/ + -JiMR-" 



(6) 



We can then calculate all the three moments of inertia A, B and 
C. 

2.2. Nonhydrostatic shape 

We here use th e observed shap e (a = 207.8 km, b = 196.7 km, 
c = 190.6 km,|^ 

omasi (1201 01) ^ to calculate the moments of in- 
ertia of Mimas, in assuming that the shape of the core is propor- 
tional to the one of Mimas, i.e. we assume 



a 



b 



c 



R ' 



(7) 



where a^., be and are the dimensions of the core, and R^ its 
mean radius (Eq[T]i. 

A quadrature over the volume of respectively the core and 
the shell gives 



Cc - III pc(^x^ + y^) dx dy dz 

^•Ju core 

= -^nacbcCc (a^ + bfjpc 

= ^jrabc{a^+b^)p,(^^\ (8) 
C. = p,{x^ +y^)dxdydz 

J J J Mimas 

- Ill Ps{^^ +y^^dxdydz 

J*J*J core 

= ^7rabc{a^+b^)p, 1 - (|^)' • (9) 

We then get C - Cc + C,,. The other moments of inertia A and B 
being obtained similarly, we have 



A = 



B = 



C = 



■^jrabc (fl^ + C-) (p,. " P-') (^) + P* 
—nabc (a^ + b^) [pc ~ P?) ("^] + P-* 



(10) 

(11) 
(12) 



We finally see that the ratio of the moments of inertia A/C - 
[b- + C-) I (fl2 + ^2) and BIC = [a^ + c^) / [a^ + b-) are inde- 
pendent on the mean radius Rc and density pc of the core, so 
every model of the internal structure of Mimas based on its ob- 
served shape (in neglecting the uncertainties on the radii a, b and 
c) will present the same rotational response. 

The interior models considered in the present study are gath- 
ered in TabU] 

3. Computing the rotation of IVIimas 

In this Section, Mimas is assumed to be a two-layers rigid body 
and the tidal contributions will be investigated in the Section |5] 



Its rotation is highly constrained by the gravitational perturba- 
tion of Saturn, and so depends on the variations of the distance 
Mimas-Saturn. That is the reason why we must understand the 
orbital dynamics of Mimas before investigating its rotation. 



3. 1 . The orbital dynamics of Mimas 

Mimas is the smallest of the main Saturnian satellites, and also 
the closest to its parent planet and t he ring s. Disc overed by 
Herschel in 1789, it is known since IStruvd ([1891 to be in 
2:1 mean-motion resonance with Tethys. More precisely, these 
two bodies are locked in an inclination-type resonance whose 
argument is 2Ai - 4At, + f2i + ^^3. the subscript 1 standing 
for the satellite S-1 Mimas, 3 for S-3 Tethys, /I, being the 
mean longitudes, and f2/ the longitudes of the ascending nodes. 
This resonance tends to raise the inclinatio ns of the sate llites 
to K 1.5° for Mimas and a; 1° for Tethys (lAllanL 1 19691) . and 
stimulates librations of the resonant argument around with 
an amplitude of x 95° and a period of w 70 years. The trap- 
ping of the system into this resonance can be explained in con- 
sidering a non-null eccentricity for Tethys that induces sec- 
ondary resonances that strong ly enhances the capture probability 
(IChampenois & Vienn^ll999 a,b). 

It is convenient to work on a Fourier-type representation of 
the orbital motion of Mimas that allows to identify every proper 
mode of the motion. The basic idea is that the variables describ- 
ing the orbital motion of Mimas can be represented as quasi- 
periodic series (and a slope for precessing angles like the as- 
cending node, the pericenter and the mean longitude), i.e. in- 
finite but converging sums of trigonometric series. The argu- 
ments of these series can be expressed as integer combination 
of a few proper modes of constant freque ncies. The existence of 
these modes comes both from the KAM (* Arnold 19631: iMoserl 



[1961 and the Nekhoroshev theories jNekhoroshevg 



1977lll979h . 



The KAM theory states that for a quasi-integrable Hamiltonian 
system (i.e. like 'K = "Hq + e'Ki where "Ho is an integrable 
Hamiltonian and eTYi a small perturbation) verifying classical 
assumptions, the motion can be considered to be on invariant 
tori (i.e. with constant amplitudes and angles depending linearly 
on time) in action-angle coordinates. For a bigger perturbation 
the Nekhoroshev theory says that the invariant tori survive over 
a timescale that is exponentially long with respect to the in- 
vert of the amplitude of the perturbation e, provided that the 
Hamiltonian of the system presents a property of steepness, that 
is an extension of the convexity. 

Such a re presentation is given by TASS1.6 ephemerides 
(IVienne & Duri ez, 1995) where the orbital motion of Mimas 
can be described using the 5 proper modes A, o), (p, ^ and 
O. A is the linear part of Mimas' mean longitude, a> is the 
main oscillation mode of the librations of the resonant argumen t 
2Ai - 4/13 +0.1+ ^li, ^ (called pi in (IVienne & Duriej. ll995h ) 
is the mean slope of A[ - 2A},, and <p - ^ and O - ^ are the mean 
slopes of respectively the longitudes of the pericenter of Mimas 
and its ascending node. The values of the frequency associated 
are gathered in Tabl3] 



3.2. Rotational model 

As for most of the natural satellites of the Solar System, Mimas 
is expected to follow the 3 Cassini Laws, or iginally described 
for the Moon (ICassiniL [1691 IColombol [T966h . i.e.: 
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Table 2. The interior models considered in the present study. The first 22 cases have been elaborated in considering Mimas as in 
hydrostatic equilibrium, while the 23'^'' is based on the observed shape. We give only one possibility for the shape model because 
the ratios of the moments of inertia A/C and B/C remain constant, so the rotational response of Mimas is the same for any interior 
model based on the shape. 



N 


Pc 


P.i 


kf 


MOI 


Ji (10--) 


C22 (10--') 


CI(mR') 


1 


1200 


800 


1.40473 


0.389636 


2.17051 


6.51152 


0.404106 


2 


1500 


800 


1.11293 


0.354953 


1.71963 


5.15889 


0.366418 


3 


2000 


800 


0.94032 


0.331801 


1.45293 


4.35878 


0.341487 


4 


2500 


800 


0.86349 


0.320705 


1.33422 


4.00267 


0.329600 


5 


3000 


800 


0.81885 


0.314001 


1.26524 


3.79571 


0.322436 


6 


3500 


800 


0.78921 


0.309439 


1.21944 


3.65831 


0.317569 


7 


4000 


800 


0.76788 


0.306100 


1.18649 


3.55946 


0.314010 


8 


1200 


1000 


1.41613 


0.390899 


2.18812 


6.56437 


0.405486 


9 


1500 


1000 


1.24455 


0.371206 


1.92301 


5.76902 


0.384026 


10 


2000 


1000 


1.17336 


0.362551 


1.81300 


5.43901 


0.374638 


11 


2500 


1000 


1.14536 


0.359061 


1.76975 


5.30925 


0.370860 


12 


3000 


1000 


1.12980 


0.357099 


1.74570 


5.23711 


0.368737 


13 


3500 


1000 


1.11969 


0.355816 


1.73009 


5.19026 


0.367350 


14 


4000 


1000 


1.11251 


0.354901 


1.71899 


5.15698 


0.366361 


15 


1200 


1100 


1.44040 


0.393565 


2.22563 


6.67688 


0.408403 


16 


1500 


1100 


1.38066 


0.386951 


2.13331 


6.39993 


0.401173 


17 


2000 


1100 


1.36451 


0.385134 


2.10836 


6.32508 


0.399189 


18 


2500 


1100 


1.35879 


0.384487 


2.09953 


6.29858 


0.398484 


19 


3000 


1100 


1.35572 


0.384139 


2.09478 


6.28435 


0.398105 


20 


3500 


1100 


1.35376 


0.383917 


2.09176 


6.27527 


0.397862 


21 


4000 


1100 


1.35239 


0.383761 


2.08963 


6.26890 


0.397692 


22 


1150.03 


1150.03 


1.5 


0.400000 


2.30951 


6.92854 


0.415397 


23 


1200 


800 


1.40473 


0.389636 


2.28639 


5.57013 


0.406273 



Table 3. The prope r frequen cies of Mimas' orbital motion (from 
TASS1.6 (IVienne&Du rieA ll995h ). 





Frequency (rad/y) 


Period (d) 


Period (y) 


A 


2435.14429644 


0.942421949 


2.580211 X 10-' 


CO 


0.08904538 


25772.62777 


70.561609 




10.19765304 


225.0452555 


0.616140 




3.81643833 


601.3285779 


1.646348 




-2.55544336 


898.0568575 


2.458746 



1 . The Moon rotates uniformly about its polar axis with a ro- 
tational period equal to the mean sidereal period of its orbit 
about the Earth. 

2. The inclination of the Moon's equator to the ecliptic is a con- 
stant angle (approximately 1.5°). 

3. The ascending node of the lunar orbit on the ecliptic coin- 
cides with the descending node of the lunar equator on the 
ecliptic. This law could also be expressed as: the spin axis 
and the normals to the ecliptic and orbit plane remain copla- 
nar 

In the case of natural satellites, they can be rephrased this way: 
the rotation of the satellite is synchronous, its angular momen- 
tum has a nearly constant inclination on an inertial reference 
plane, and is located in the plane defined by the normal to the 
orbital plane and to the Laplace Plane. The Laplace Plane is the 
plane normal to the rotation axis of the orbital frame, i.e. it is 
defined with respect to the orbital precessional motion. It has the 
property to minimize the variations of the orbital inclinations. 
For satellites orbiting close to their planet as it is the case here, 
the equatorial plane of Saturn is so close to the Laplace Plane 
that it can be used for describing the rotational dynamics. 

Our rotational model is similar to the one already used in 
e.g. (iNovelles et al.Ll2008tlNoveUesll2010h for studying the rigid 
rotation of the Saturnian satellites Titan, Janus and Epimetheus. 




Fig. 2. The An doyer variables (partially reproduced from 
(lHenrardLl2005ah '). 



We consider Mimas as a rigid triaxial body whose matrices 
of inertia reads 



A ^ 

B 

00c 



(13) 



with A<B <C. 

The dynamical model is a 3 -degree of freedom one. We will 
use the Andoyer variables which requires a decomposition with 
3 references frames : 



1. 



An inertial reference frame (ci, 62,63). We used the one in 

which the orbital ephemerides are given, i.e. mean Saturnian 

equator and mean equinox for J2000.0 epoch. 

A frame («i,«2,«3) bound to the angular momentum of 

Mimas. 

A frame (/i,/2,/3) rigidly linked to Mimas. 



W e first use Andoyer's variables dAndoved, 119261: iDepritl 
II967h . which are based on two linked sets of Euler's angles. The 
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first set (h, K, g) locates the position of the angular momentum in 
the first frame {61,62,63), while the second one, (g, J, I), locates 
the body frame (/i, /z, /s) in the second frame tied to the angular 
momentum (see Fig.|2|. 

The canonical set of Andoyer's variables consists of the three 
angular variables l,g,h and their conjugated momenta L,G,H 
defined by the norm G of the angular momentum and two of its 
projections: 



L — G cos J, 
G, 

H = GcosK. 



(14) 



Unfortunately, these variables present two singularities: 
when / = (i.e., the angular momentum is colinear to /j), / 
and g are undefined, and when K - (i.e., when Mimas' prin- 
cipal axis of inertia is perpendicular to its orbital plane), h and g 
are undefined. That is the reason why we shall use the modified 
Andoyer's variables: 



p - I + g + h, 
r — —h. 



^9 = V^sin^, 



= = p(i - cos/:), 

= 2P sin^ f , 



^1 = V^cos^, 



(15) 

where n is the body's mean orbital motion , q — -I, and Q - 
G - L - G{ \ - cos J) = 2G sin^ j. With these new variables, the 
singularity on I has been dropped. Using these variables has a 
great mathematical interest, because they are canonical, so they 
simplify an analytical study of the system, as was done in the 
previous works mentioned above. Our study here is essentially 
numerical, but we keep these variables, in order to be consistent 
with previous studies. We later derive other output variables, that 
are more relevant from a physical point of view. 

In these variables, the kinetic energy T - jCo-G of the system 
reads: 



ri + 72 



ri - 72 1-71+72 



71 - 72 21 



(16) 



with 



with 



and 



V, 



pi 



-2C^[7i(-, 



■riixl-yl)] 



(20) 



Vp2 = - 



15 gM„^R„ 

-jCJ2p 



73 



X [71(4 + 72(4 -3'p)]' 



(21) 



where is the gravitational constant, Mp the mass of the per- 
turber, J2p its J2, Rp its mean radius, dp the distance between 
the perturber's and Mimas' centers of mass, and Xp and yp the 
two first components of the unit vector pointing to the center of 
mass of the perturber, from the center of mass of the body, in the 
reference frame {f\,f2,fi)- Vp\ expresses the perturbation due 
to a pointmass perturber, while Vp2 represents the perturbation 
due to its J2, assuming that the body is in the equatorial plane 
of the perturber. As shown in (Henrard, 2005c), it is a good ap- 
proximation if the sine of the angle between Saturn's equatorial 
plane and the orbit is small. In the case of Mimas, this angle (i.e. 
Mimas' orbital inclination) is ^ 1.5° ~ 2.6 x 10"^ rad, so we 
can consider that its sine is always smaller than 3 x 10"^. This 
assertion also assumes that the obliquity of Mimas is very small, 
what will be checked in this study. 

Usually the orbital ephemerides give us the location of the 
perturber in the inertial frame, so we have to perform 5 rotations 
to convert the coordinates from the inertial frame to (/i,/2,/3). 
More precisely, if we name (.v,,^,, z,)^ the unit vector locating 
the perturber's center of mass in the inertial frame, we have 



Xp ^ 

yp 

y Zp ) 



■■ R3(-l)Ri{-J)Ri{-g)Ri(-K)R,{-h) 



y. 

\ Zi j 



(22) 



with 



^3(0) 



and 



R^ 



cos - sin ' 
sin <p cos <p 
1, 



no ^ 

COS ^ - sin ( 
sin0 cos(i 



(23) 



(24) 



and 



IC-A-B MR^ 
71 = — = Ji- 



2C 



C 



72 



B-A MR^ 



2C 



C 



(17) 



(18) 



In these last 3 formulae, a> is the instantaneous vector of rota- 
tion, M is the mass of Mimas, R its mean radius, and J2 and C22 
the two classical normalized gravitational coefficients related re- 
spectively to the oblateness and equatorial ellipticity of the con- 
sidered body. 

The gravitational disturbing potential due to an oblate per- 
turber p reads (.Henrard. 2005c): 



Vp = y„i + Vp2 



(19) 



Finally, the total Hamiltonian of the problem reads: 



H 



nP n r o On 



71 +72 



-7: ' 
3&M^,. 5 R 



71 -72 



2« di 



71 -72""' 1 -71 +72 



x[7i^ +y]^) + j2{x\^-y\^)]. 



(25) 



where the index stands for Saturn. We will use this 
Hamiltonian for a numerical study of the rotation. An analyti- 
cal study can show that the Hamiltonian ( |25] | can be reduced to 
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'Hiu, V, w, U, V, W) ^u„U + (jjyV + w„W + V{u, v, w, U, V, W) 

(26) 

where !P represents a perturbation, and the three constants co,,,, 
(jL>y and (JL>„ are the periods of the free oscillations around the 
equilibrium defined by the Cassini Laws. This last Hamiltonian 
is obtained after several canonical transformations, the first one 
consisting in expressing the resonant arguments cr - p - A + n 
and p = r + f2 respectively associated with the 1 : 1 spin-orbit 
resonance and with the orientation of the angular momentum, A 
and f2 being the orbital variables defined above. The complete 
calculatio n is beyond the s co pe of this paper, the re ader can find 
details in (lHenrardil2005allbtlNovelles et al.Ll2008h . 

3.3. A numerical study 

In order to integrate numerically the system, we first express 
the coordinates of the perturber {xf^,y-^) with the numerical 
ephemerides and the rotations given in (Eql22ll. in the body 
frame (/i,/2,/3). As explained bef ore, the ephemerides ar e 
given by the TASS1.6 ephemerides dVienne & Duriezi Il995l) . 
This way, we get coordinates depending of the canonical 
variables. Then we derive the equations coming from the 
Hamiltonian dZSl i: 



dp 


dH 


dP 


dH 


It ' 


" dP' 


dt 


dp' 


dr 


dH 


dR 


dH 


dt ' 


' dR' 


'dt " 


dr' 


d^g _ 


dH 


driq 


dH 


dt 




dt 





(27) 



We integrated over 200 years using the Adams-Bashforth- 
Moulton 10th order predictor-corrector integrator The solutions 
consist of two parts, the forced one, directly due to the perturba- 
tion, and the free one, that depends on the initial conditions. The 
initial conditions should be as close as possible to the exact equi- 
librium, that is assumed to be the Cassini State 1 in 1 : 1 spin-orbit 
resonance, to have low amplitudes of the free libr ations. For that, 
we h ave used the iterative algorithm NAFFO dNovelles et al.L 
l20Tlh to remove the free librations from the initial conditions, 
after they have been identified by frequency analysis. 

The frequency analysis algorithm we used is based on 
Laskar's original idea, named NAFF as Num erical Analysi s of 
the Fundamental Frequencies (s ee for instance lLaskarl (Il993l) for 
the method, and iLaskad (l2005l) for the convergence proofs). It 
aims at identifying the coefficients at and ojk of a complex sig- 
nal /(f) obtained numerically over a finite time span [—T\ J] and 
verifying 

n 

fit) ^ exp((WA.f), (28) 

A-=l 

where wj- are real frequencies and complex coefficients. If the 
signal /(f) is real, its frequency spectrum is symmetric and the 
complex amplitudes associated with the frequencies and -wa 
are complex conjugates. The frequencies and amplitudes associ- 
ated are found with an iterative scheme. To determine the first 
frequency coi, one searches for the maximum of the amplitude 
of 

<^(w) /(f), exp(;wf) >, (29) 



where the scalar product < /(f), g(t) > is defined by 

< /(f), g(t) >= ^ f(t)^x(t)dt, (30) 
and where ;t'(f) is a weight function, i.e. a positive function with 

^j^X(t)dt^l. (31) 

Once the first periodic term exp((Wif) is found, its complex am- 
plitude fli is obtained by orthogonal projection, and the process 
is started again on the remainder / (f) = /(f) - oi exp{icoit). The 
algorithm stops when two detected frequencies are too close to 
each other, which alters their determinations, or when the num- 
ber of detected terms reaches a maximum set by the user. This 
algorithm is very efficient, except when two frequencies are too 
close to each other In that case, the algorithm is not confident 
in its accuracy and stops. When the difference between two fre- 
quencies is larger than twice the frequency associated with the 
length of the total time interval, the determination of each fun- 
damental frequency is not perturbed by the other ones. Although 
the iterative method suggested by Champeno is (1998) allows to 
reduce this distance, some difficulties remain when the frequen- 
cies are too close to each other. 

3.4. Outputs 

In order to deliver theories of rotation that can be easily com- 
pared with observations, we chose to express our results in the 
following variables; 

- Longitudinal librations, 

- Latitudinal librations, 

- Orbital obliquity e (the orientation of the angular momentum 
of Mimas with respect to the normal to the instantaneous 
orbital plane), 

- Motion of the rotation axis about the pole axis. 

There are at least two ways to define the longitudinal libra- 
tions. We can for instance consider the librations about the exact 
synchronous rotation, i.e. p- < n > f, usually called physical 
librations . In this case, we have used for < n > the frequency 
associated with the proper mode A, i.e. Mimas' mean longi- 
tude. Another way to consider the longitudinal librations is to 
work on the librations about the Mimas-Saturn direction. We 
will call these librations tidal librations because they represent 
the misalignment of the tidal bulge of the satellite. The differ- 
ence between these two librations is known as optical librations , 
only due to the orbital motion of Mimas around Saturn. The 
reader can find graphica l descriptions of these librations in 
iMurrav & Dermottid 19991) . Fig.5.16. 

The latitudinal librations are the North-South librations of 
the large axis of the considered body in the saturnocentric refer- 
ence frame that follows the orbital motion of the body. They are 
analogous to the tidal librations that are the East- West librations. 
In order to get the tidal longitudinal librations and the latitudi- 
nal librations, we first should express the unit vector /i (i.e. the 
direction of Mimas' long axis) in the inertial frame {ei,e2,e^). 
From (Eql22li and the definitions of the Andoyer modified vari- 
ables (EqfTSll. we get: 
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/l 



- (cos r(cos(/7 + r - I) cos / - sm{p + r - I) cos J sin /) 

+ sin r(cos K(sm(p + r - I) cos / 

+ cos(/7 + r — 1) cos y sin /) - sin K sin / sin I))ei 
+ (- sin r{cos{p + r — I) cos Z - sin(/7 + r - Z) cos 7 sin /) 

+ cos r(cos K{s,\n{p + r - I) cos / 

+ cos(;9 + r - I) cos J sin /) - sin K sin J sin 0)^2 
+ (sin K{sm{p + r - I) cos / + cos(p + r - I) cos 7 sin Z) 

+ cos /iT sin J sin 0^3 ■ (32) 



The tidal longitudinal librations iff and the latitudinal ones are 
found this way: 



and 



nfi 



(33) 



(34) 



where n is the unit vector normal to the orbit plane, and t the 
tangent to the trajectory. We get these last two vectors by: 



and 



XXV 

\\x X vll 



nxx 

\\n X jcll 



(35) 



(36) 



where x is the position vector of the body, and v its velocity. 

Finally, the motion of the rotation axis about the pole is de- 
rived from the wobble /, it is given by the two variables Qi and 
Q2 defined as: 



and 



Qi - sin J sin l^l + 
Q2 = sin / cos l^l + 



■/2 + 2C22 ' 
C 



c 



(37) 



(38) 



they are the first two components of the unit vector pointing 
at the instantaneous North Pole of Mimas' rotation axis, in the 
body frame of Mimas. These quantities are fi nally mul t iplied by 
the polar radius of the satellite, i.e. 190.6 km (lThomasll2010l) to 
get a deviation in meters. 



4. Results 

We here present the outputs of our numerical study of the rota- 
tion of Mimas. We first give the example of a non-hydrostatic 
model of Mimas based on its observed shape, then we compare 
the results with the rotational response of the first 22 models of 
Tab|2] obtained in considering Mimas to be in hydrostatic equi- 
librium. 



4.1. Non-hydrostatic Mimas based on its sinape 

As already mentioned, this case is unique, because changes in 
the size of the core do not affect the ratios of the moments of 
inertia A/C, and BjC, and the coefficients yi and 72 (EqfTTland 
[TS] ). As a consequence, there is a unique rotational behavior of 



Table 4. Frequencies and periods of the free librations of Mimas, 
in the shape model. These values have been obtained numeri- 
cally. 



Proper 


Frequency 


Period T 


mode 


(rad/d) 


(d) 


u 


2.704622 


2.323129 


V 


0.778015 


8.075914 


w 


0.621287 


10.113182 



Mimas for any homogenous or 2-layer model using this specific 
model based on the observed shape. 

The free librations around the equilibrium are assumed to be 
damped, it is anyway important to know their frequencies w„, w,, 
and u„ (or periods r„, T,, and T,,,) because they characterize the 
way the system reacts to external sinusoidal excitations, that are 
here due to the variations of the distance between the Sun and 
Mimas. 

The frequencies of the free librations are listed in TabHl The 
proper mode u roughly represents the free longitudinal libra- 
tions, V the free librations of the obliquity, and w the wobble, i.e. 
the free polar motion of Mimas. These frequencies have been 
deduced from the frequency analysis of the modified Andoyer 
variables (cf.EqfTSll. 

So, the proper modes involved in the Fourier representations 
of the librations of Mimas are the forced modes due to the or- 
bital motion of Mimas around Saturn (cf.TabO and the free 
ones (Tab|4|i. The arguments of the sinusoidal components of 
the quasi-periodic decompositions of the variables of the rota- 
tion are integer combinations of these proper modes. If we con- 
sider that the free librations are damped, the solutions should be 
only composed of the forced modes. 

The forced librations of Mimas modeled from its observed 
shape are given in Tab|5]to|8] These tables give the solutions un- 
der the form of periodic time series, in cosines. We can see that 
the main difference between the physical and the tidal librations 
is in the presence in the physical librations of a long-period term 
(rs 70 years) with a high amplitude (« 43°, i.e. ^ 86° peak-to- 
peak) due to the librations of the argument of the orbital reso- 
nance between Mimas and Tethys. As explained in Rambaux et 
al. (2010) and (201 1), the amplitude of the long period librations 
are equal to the magnitude of the orbital perturbations because 
at long period the body is oriented toward the central planet. As 
a consequence, by analysing the tidal librations, the long period 
librations vanish. There is also a large difference in the ampli- 
tude given for the tidal and physical longitudinal librations. As 
explained above, this difference is due to optical librations, with 
amplitude 2e ^ 3.8 x 10"^ rad ^ 2.2°. 

The latitudinal librations of Mimas (Tab|7]i are significantly 
smaller 2 arcmin vs. 2.5° for the tidal longitudinal librations), 
and so could hardly be used in the framework of observations of 
the rotation of Mimas (except if there are free oscillations due to 
a recent unexpected excitation). The mean obliquity of Mimas 
(Tabl8]l is of the same order of magnitude. 

4.2. For a tiydrostatic Mimas 

We performed the same numerical study of the 22 hydrostatic 
configurations of Mimas given in Tab|2] The results are gathered 
in Tab|9l In this table, the amplitudes of the tidal and physical 
longitudinal librations indicated are related to the mode A-(p + ^ 
(period: 0.944898 d), while the latitudinal ones are related to the 
mode A + (p - ( (period: 0.939962 d). The main physical reason 



8 B. Noyelles et al.: The rotation of Mimas 

Table 5. Forced tidal longitudinal librations of Mimas, in the shape model. The series are in cosine. 











Frequency 


Period 


Amplitude 


Phase 


A 




4> 


( 


(rad/y) 


(d) 


(arcmin) 


at J2000 


1 




-1 


1 


2428.763080 


0.944898 


157.73363 


-79.177° 


1 


1 


-1 


1 


2428.852395 


0.944863 


5.72739 


-116.201° 


1 


-1 


-1 


1 


2428.673643 


0.944933 


4.05163 


137.568° 


2 




-2 


2 


4857.526150 


0.472449 


1.83667 


-68.313° 


1 






1 


2438.960801 


0.940947 


1.32391 


-148.065° 



Table 6. Forced physical longitudinal librations of Mimas, in the shape model. The series are in cosine. 









Frequency 


Period 


Amplitude 


Phase 


A w 




-r 


(rad/y) 


(d) 


(arcmin) 


at J2000 


1 






0.08904538 


25772.62777 


43.61° 


51.354° 


- 3 






0.26713614 


8590.87592 


43.261 arcmin 


-25.913° 


1 


-1 


1 


2428.763080 


0.944898 


26.075 arcmin 


101.355° 




1 




10.19765304 


225.04526 


7.828 arcmin 


-157.744° 


1 


-1 




10.10860766 


227.02728 


3.657 arcmin 


-119.032° 


1 


1 




10.28669842 


223.09718 


3.532 arcmin 


-16.309° 



Table 7. Forced latitudinal librations of Mimas, in the shape model. The series are in cosine. 











Frequency 


Period 


Amplitude 


Phase 


A 


a) 




( 


(rad/y) 


(d) 


(arcmin) 


at J2000 


1 




1 


-1 


2441.516177 


0.939962 


2.07096 


77.130° 


1 


1 


1 


-1 


2441.605507 


0.939928 


0.06829 


39.984° 


1 


-1 


1 


-1 


2441.426665 


0.939997 


0.0414 


-66.603° 



of these librations is the variations of the distance Mimas-Saturn 
during an orbital period. 

In order to make the results more readable, we present them 
graphically in Figl3] The plots present a clear dependency of 
the amplitudes of librations on the densities of the core and the 
shell. We can in particular notice that the longitudinal (i.e. tidal 
and physical) librations have a larger amplitude when the density 
of the core is lower, it is due to the fact that a concentration of 
the mass in the core lowers the moments of inertia of the body, 
and so tends to limit its amplitude of response to sollicitations. 
Finally we can see that the dependency on is small for - 
llQQkg/in^, it is because in this case, ps is close to the mean 
density of Mimas (i.e. 1 150.03kg /m^), as a consequence the core 
is small and Mimas is close to be homogeneous. 

4.3. A small polar motion 

As for the other outputs, we present the forced polar motion 
of Mimas (i.e. after removal of the free wobble) as a sum of a 
trigonometric series (TabfTOt. We can see that this motion is ex- 
pected to be small, the highest amplitude being a; 15 meters. The 
sum of all these amplitudes can reach 40 meters, so we can con- 
sider these 40 meters as the upper bound of the polar motion. An 
analysis of the polar motions for the different hydrostatic Mimas 
do not exhibit significant differences. 



5. Tidal Dissipation 

This section is dedicated to study the influence of the tidal torque 
on the rotational motion of Mimas. We introduce the tidal torque 
in a Lagrangian formalism and follow the approach of Williams 
et al. (2001) and used recently in Rambaux et al. (2010) and 
Robutel et al. (2011). The starting equation is the angular mo- 



mentum equation 

— +(oaG^T (39) 
at 

where a> is the angular velocity vector, the angular momentum 
G - Icj with / the tensor of inertia, and T is the external gravi- 
tational torque expressed as 

3GMu 

T = r— A lu (40) 

where u is the cosine director of Saturn in the reference frame 
tied to Mimas, and My^ its mass. 

The dissipation is due to the tidal and centrifugal potentials 
that deform the satellite. In this case, the tensor of inertia / be- 
comes a constant plus a time-variable part resulting from the de- 
formation. The time-variable part does not react instantaneously 
and therefore presents a time delay 6t characteristics of the rhe- 
ological properties of the body (see section|2]i. 

In addition, the dynamical equation of the rotational motion 
Eq. [39] may be linearized by using the synchronous spin-orbit 
resonance of the body implying that (jJi,(jl>2 « W3 ~ n and 
M2,M3 << u\ ~ \ where ui,U2,U3 are the coordinates of the 
cosine director along the principal axis of inertia of Mimas. 

By using these approximations and focusing on the libration 
in longitude, the main tidal torque is expressed as (Williams et 
al. 2001) 

3GM2 

T,^-k2R'^^(UuU;,-UnUu), (41) 

where Uij = (7) ^^^^ indicates the time delay part. 

Then, we used the same approach than in Rambaux et al. 
(2010). We introduce the rotation angle ip similar to the sum of 
the Andoyer angles / + g because the polar motion is small as 
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A 


w 





O 


Frequency 
(rad/y) 


Period 
(d) 


Amplitude 
(arcmin) 


Phase 
at J2000 













OO 


2.13468 




2 




1 


-1 


4883.032354 


0.469981 


0.07372 


-27.144° 


2 


-1 


1 


-1 


4882.943311 


0.469990 


0.07179 


14.019° 


2 


1 


1 


-1 


4883.121227 


0.469973 


0.06866 


1 19.800° 




2 






0.177998 


12893.06 


0.04817 


103.293° 



Table 9. Periods of the free librations and amplitudes (in arcmin) of the forced librations for the different models assuming that 
Mimas is at the hydrostatic equihbrium. 





T, 




T„ 


Tidal 


Latitudinal 


Mean 


Physical 


N 


(d) 


(d) 


(d) 


librations 


librations 


obliquity 


librations 


1 


2.143878 


7.885550 


11.621674 


163.398 


2.016 


2.086 


31.744 


2 


2.294081 


8.997072 


13.222674 


158.577 


2.314 


2.384 


26.914 


3 


ZAOl 1 1 1 


9.908107 


12.763086 


155.609 


2.559 


2.631 


23.944 


4 


2.468518 


10.416095 


13.618627 


154.248 


2.693 


2.765 


22.582 


5 


2.507237 


10.742236 


14.181571 


153.442 


2.780 


2.853 


21.776 


6 


2.534519 


10.975456 


14.591169 


152.900 


2.843 


2.917 


21.234 


7 


2.555063 


11.152890 


14.906138 


152.508 


2.891 


2.966 


20.841 


8 


2.138844 


7.849766 


11.569966 


163.583 


2.006 


2.076 


31.922 


9 


2.220696 


8.443608 


12.426948 


160.777 


2.166 


2.234 


29.115 


10 


2.257477 


8.730512 


12.839826 


159.593 


2.243 


2.313 


27.930 


11 


2.274983 


8.851150 


13.013162 


159.124 


2.275 


2.346 


27.461 


12 


2.284054 


8.920304 


13.112486 


158.862 


2.293 


2.363 


27.199 


13 


2.290031 


8.966015 


13.178106 


158.692 


2.306 


2.378 


27.028 


14 


2.294326 


8.998974 


13.225425 


158.571 


2.315 


2.384 


26.907 


15 


2.128297 


7.775062 


11.461869 


163.975 


1.986 


2.054 


32.314 


16 


2.154684 


7.962727 


11.733240 


163.008 


2.037 


2.104 


31.347 


17 


2.162067 


8.015705 


11.809777 


162.745 


2.051 


2.119 


31.084 


18 


2.164709 


8.034707 


11.837202 


162.652 


2.056 


2.124 


30.991 


19 


2.166135 


8.044989 


11.852091 


162.602 


2.059 


2.127 


30.941 


20 


2.167044 


8.051527 


11.861492 


162.570 


2.060 


2.132 


30.909 


21 


2.167685 


8.056158 


11.868221 


162.548 


2.062 


2.130 


30.886 


22 


2.106951 


7.625231 


11.244996 


164.792 


1.946 


2.014 


33.132 



Table 10. Polar motion of Mimas qi + iq2, in the shape model. The series are in complex exponentials. 











Frequency 


Period 


Amplitude 


Phase 


A 


a) (p 




O 


(rad/y) 


(d) 


(m) 


at J2000 


1 




-1 




2441.516177 


0.939962 


15.277 


77.130° 


-1 




1 




-2441.516177 


0.939962 


14.368 


102.870° 








1 


-12.7530964 


179.951077 


4.498 


113.654° 








-1 


12.7530964 


179.951077 


3.441 


66.208° 


1 


1 1 


-1 




2441.605507 


0.940313 


0.503 


40.052° 


2 




2 


-1 


4870.2792464 


0.471212 


0.485 


87.878° 


-1 


-1 -1 


1 




-2441.605507 


0.940313 


0.473 


139.948° 



shown in Figure |2] where J and (the nutation angle) are small. 
The libration angle y is defined as = M + y representing the 
physical libration in longitude of the body. The cosine director 
U2 ~ i - y is of the order of the difference between the orbital 
variation s and the physical libration y. We note that U2 corre- 
sponds to the tidal libration if/ defined in section 13.41 and their 
amplitude is small as shown in Table 5. The quantity s, the or- 
bital variation, is defined as the difference between the true and 
the mean longitude of Mimas and represents the oscillation of 
the orbital longitude of Mimas that may be expressed in Fourier 
series as 



s = ^Hi sin (a>it + a,). (42) 

i C Gm 



Then, by developing u\ and in Taylor series for small dt, the 
dynamical equation becomes; 

3 G ^^^h 

Cy + -{B - A)-^ sin2(y - s) ^ -k2R' -^6t{y - s). (43) 

2 a" 

As shown in the previous section, the quantity y - s is always 
small (see Table 5) allowing to simplify the sine function by its 
angle. In addition, the eccentricity of Mimas is small and so a/r 
is equal to 1 at first order in eccentricity. Finally, we obtain a 
forced dissipative harmonic oscillator with the frequency ojq = 
« V3(B - A)/C and the dissipative rate A expressed as 

3k2R^ n' 

2 A = — 6t. 44) 
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Fig. 3. Librations of hydrostatic Mimas, depending on the densities of the core and the shell. 



m being the mass of the satellite and k2 is the Love number of 
Mimas. 

As for the conservative case, the amplitude of terms associ- 
ated with the long period is almost equal to the magnitude of the 
oscillation s. The solution may be expressed as 



After a restricted analytical approach to validate the numeri- 
cal results, we discuss the possibility to observe the rotation 
of Mimas and in particular to discriminate the different interior 
models. Then we focus on the non-hydrostatic contributions. 



7 



Ad sin {codt + (pd)e''^' + ^ x-, cos (w,f + a{)+yi sin (w,f + a,), 6.1. Analytical approach 



(45) 

where Ad and (pd are constants of integration. The first term de- 
cays with time scale IjA and its resonant frequency is Ud - 

■yjoJ^ - /l^. The periodic term of the particular solution y is com- 
posed of the in-phase y,- and out-of -phase x, terms 



yi = Hi 



H, 



(46) 

At first order, the expression of x, may be simplified as 



Xi = -0.9054-^//,- 



(47) 



expressed in radians and by introducing the dissipation factor 
as St = (nQT^. For short period librations at 0.944898 days 
the Xi is 1.32 mas with ^ = 10"^ (this is the value used by 
iMever & WisdomI (l2008h ) and the resulting displacement at the 
surface of the satellite at the periaster passage is also negligible 
0.0013 m. The damping time 1 /A is about 6,000 years. If we con- 
sider ^ to be 100 times bigger, i.e. closer to the expected value 
of Enceladus, we have a displacement at the surface of ^ 0.13 m 
and a damping time of x 60 years. For librations at long period 
Xj is definitely negligible because w, is very small. 

6. Discussion 

One of the aims of this theoretical study is to prepare the in- 
terpretation of potential observations of the rotation of Mimas. 



We here compare with classical analytical formulae for the main 
term of the physical and tidal longitudinal librations and the 
mean obliquity, for which deriving accurately these amplitudes 
is quite straightforward. 

6.1.1. Longitudinal librations 

An analytical study of the longitudinal librations of a satellite in 
1:1 spin-orbit resonance and on a kepler ian orbit can be found 
for instance in iMurrav & DermottI d 19991) . Let us call if/ the am- 
plitude of the main term (i.e. associated with the mode A + (p-0 
of the tidal librations, and y for the physical ones. We have from 
,Murrav & Dermotti(il999i) : 



if, = 



r = 



and 



-2e 



2e 



3 1 -5e^ + —e*\^ 12 ^ 

C \ 16 CfimR^) 



(48) 

(49) 
(50) 



1 - 5e + — e 
16 

(51) 



e being the eccentricity of Mimas. We can see that this amplitude 
is bigger when the ratio is closer to unity, or when 1 ICni is closer 
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to C/{mR^). We can see from the Tab|2]that C22 is of the order 
5 X 10"^ while C ^ QAmR^, i.e. C22/C/{mR^) ^ 1/80. Thus, the 
ratio i^ c/(mR-) closer to unity when C22 is bigger, what is the 
case for the smallest values of pc- The Fig|3]confirms this trend, 
while the TabfTTI settles the validity of the analytical formulae 
(g8]i and (gUi. 



6.1.2. Mean obliquity 

We here use the ana lytical study of IWard & Hamiltonl (l2004t) 
(see iNovellesI (1201 Ol) for an application to natural satellites in 
spin-orbit resonances) for the location of the Cassini States. 
Mimas is expected to be locked at the Cassini State 1, i.e. the 
most stable one, characterized by; 



sin/ 



3,1 J2+2C22 

2^ C/(mR2) 



+ COS / 



(52) 



e being the mean obliquity of Mimas, Cl the precessional rate of 
its orbital ascending node, and / its inclination on the Laplace 
Plane, here assumed to be the equator of Saturn at J2000. 

From the definition the orbital proper modes of Mimas, 
we can approximate f2 by ~ ^ = -6.37188169 rad/y, this 
yields a regressional period of 360.1657 days. In assuming J2 ~ 
2 X IQ-^, C22 ~ 6 X 10-^ and C ^ OAmR^ from Tabl2l we 
have 4^7777^^ ~ -37.26 while sin / is very small and cos/ 

close to unity (the mean inclination of Mimas / being of the or- 
der of 1.5° = 2.6 X 10"^ rad). So, bigger values of the quantity 
J2 + 2C22 will yield a smaller obliquity. Once more, these values 
are reached for the lowest estimations of p^, the Figl3] confirm- 
ing this tendency. The validity of the analytical formula ( l52l i is 
checked in TabfTTI 

The analytical validation meets the following trouble: how 
to evaluate the mean eccentricity and inclination required in the 
analytical formulae, i.e. how to average them? These formulae 
have been derived in assuming a Keplerian orbit, while the or- 
bit of Mimas is perturbed by the oblateness of Saturn and the 
mutual interactions with the other satellites, inducing an orbital 
resonance with Tethys. As a consequence, its eccentricity and 
inclination are far from constant. 

We have, from lVienne & Duriei (Il995h : 



z(t) = e(t)exp(im(t)) 

= 1.59817 X 10^2 exp(((6.38121472f+ 356.521°)) 

+ 7.2147 X 10"^^ exp {i (6.29216934f + 137.197°)) (53) 

+ 7.11 14 X 10"^ exp (i(6.47026010f+ 35.846°)) -H..., 

at) = sinj^j(;f2(0) 

= 1.18896 X 10^2 exp (((-6.37188169f + 234.213°)) 

+ 5.3177 X 10"^^ exp (( (-6.46092707f + 14.888°)) (54) 

+ 5.3017 X 10"^^ exp (((-6.28283631fH- 273.538°)) -H..., 



the frequencies being in rad/year, and the time origin J 1980. As 
we can see, the mean eccentricity should be at least a: 1 .6 x 10"^, 
probably higher (same for the mean inclination, that should be 
at least ^ 1.4°). In the TabUI] we use e = 1.92 x 10"^ and 
/ - 1.68°, this arbitrary choice minimizes the relative errors and 
is consistent with the TASS1.6 theory. 



6.2. Observational possibilities 

It would be challenging to constrain the orientation and interior 
structure of Mimas using its rotation. The first expected result 
is the confirmation that Mimas is in the Cassini State 1 with the 
1:1 spin-orbit resonance. Another challenge would be to detect 
the longitu dinal hbration s, that have been actually obser ved for 
the M oon jKoziell (Il967h ). the Martian satellite Phobos dBurns 
(1972)), and the Saturnian satellite Epimetheus dTiscareno et al. 
(l2009h V To estimate the required accuracy of the observations, 



we convert the rotation outputs into kilometres (TabfTSli. 



Table 12. Expected librations and mean obliquity of Mimas, in 
km. The mean obliquity e has been multiplied by the polar radius 
c = 190.6 km, while the librations have been multiplied by the 
Saturn-facing radius a = 207.8 km. The case 23 is derived from 
the shape model. 





Physical 


Tidal 


Latitudinal 


Mean 


N 


librations 


librations 


librations 


Obliquity 


1 


1.919 


9.877 


0.122 


0.116 


2 


1.627 


9.585 


0.140 


0.132 


3 


1.447 


9.406 


0.155 


0.146 


4 


1.365 


9.324 


0.163 


0.153 


5 


1.316 


9.275 


0.168 


0.158 


6 


1.284 


9.242 


0.172 


0.162 


7 


1.260 


9.219 


0.175 


0.164 


8 


1.930 


9.888 


0.121 


0.115 


9 


1.760 


9.718 


0.131 


0.124 


10 


1.688 


9.647 


0.136 


0.128 


11 


1.660 


9.618 


0.138 


0.130 


12 


1.644 


9.603 


0.139 


0.131 


13 


1.634 


9.592 


0.139 


0.132 


14 


1.626 


9.585 


0.140 


0.132 


15 


1.953 


9.912 


0.120 


0.114 


16 


1.895 


9.853 


0.123 


0.117 


17 


1.879 


9.837 


0.124 


0.117 


18 


1.873 


9.832 


0.124 


0.118 


19 


1.870 


9.829 


0.124 


0.118 


20 


1.868 


9.827 


0.125 


0.118 


21 


1.867 


9.825 


0.125 


0.118 


22 


2.003 


9.961 


0.118 


0.112 


23 


1.576 


9.534 


0.125 


0.118 



As expected, the longitudinal librations are significantly big- 
ger (a few kilometres) than the mean obliquity and the latitudinal 
librations (with an amplitude smaller than 200 meters). The am- 
plitude of the librations given are related to the quasi-periodic 
decompositions, so the peak-to-peak amplitudes are twice big- 
ger The reader should keep in mind that the physical and tidal 
librations are two expressions of the same quantity, so are not in- 
dependent. We can consider that the detection of the longitudinal 
librations would require an accuracy of about 1 km, while using 
them to invert the internal structure of Mimas would require an 
accuracy at least ten times better 

6.3. Non-hydrostatic contributions 

The study of the non-hydrostatic Mimas, based on the shape 
model, does not exhibit a significant possibility to discriminate 
a non-hydrostatic Mimas from a hydrostatic one from observa- 
tions. This is not surprising considering Mimas' nearly hydro- 
static global shape. But a non-hydrostatic Mimas could result in 
an offset between the ellipsoid o f shape and the ellips oid of iner- 
tia, as investigated for Janus bv lRobutele"tan(l201lh . for which 
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Table 11. Analytical confirmation of the numerical results given in Tab|2] The analytical formulae used are Eq|5Tl|49]|48]&|52] the 
obtained values being compared with the ones given in Tab|9l 





Tu 




r 


Ay 






6 


Ae 


N 


(d) 




(arcmin) 




(arcmin) 




(arcmin) 




1 


2.145171 


0.060% 


31.572 


0.542% 


163.581 


0.112% 


2.088 


0.104% 


2 


2.294910 


0.036% 


26.778 


0.506% 


158.787 


0.133% 


2.397 


0.547% 


3 


2.410238 


0.102% 


23.825 


0.497% 


155.835 


0.145% 


2.651 


0.741% 


4 


2.471008 


0.101% 


22.471 


0.493% 


154.480 


0.150% 


2.790 


0.890% 


5 


2.509751 


0.100% 


21.669 


0.490% 


153.679 


0.154% 


2.880 


0.956% 


6 


2.537080 


0.101% 


21.131 


0.487% 


153.140 


0.157% 


2.945 


0.966% 


7 


2.557614 


0.100% 


20.740 


0.487% 


152.749 


0.158% 


2.994 


0.960% 


8 


2.140163 


0.062% 


31.755 


0.522% 


163.765 


0.111% 


2.078 


0.108% 


9 


2.221696 


0.045% 


28.965 


0.514% 


160.975 


0.123% 


2.243 


0.411% 


10 


2.259963 


0.110% 


27.788 


0.508% 


159.798 


0.128% 


2.323 


0.429% 


11 


2.275851 


0.038% 


27.321 


0.508% 


159.331 


0.130% 


2.356 


0.446% 


12 


2.284904 


0.037% 


27.061 


0.507% 


159.071 


0.131% 


2.376 


0.537% 


13 


2.290872 


0.037% 


26.892 


0.505% 


158.901 


0.132% 


2.388 


0.438% 


14 


2.295156 


0.036% 


26.771 


0.506% 


158.780 


0.132% 


2.398 


0.569% 


15 


2.129674 


0.065% 


32.145 


0.522% 


164.155 


0.110% 


2.057 


0.170% 


16 


2.155925 


0.058% 


31.184 


0.522% 


163.193 


0.113% 


2.110 


0.266% 


17 


2.163275 


0.056% 


30.922 


0.520% 


162.932 


0.115% 


2.124 


0.251% 


18 


2.165905 


0.055% 


30.830 


0.520% 


162.839 


0.115% 


2.130 


0.263% 


19 


2.167325 


0.055% 


30.780 


0.520% 


162.789 


0.115% 


2.132 


0.256% 


20 


2.168230 


0.055% 


30.748 


0.520% 


162.758 


0.116% 


2.134 


0.106% 


21 


2.168868 


0.055% 


30.726 


0.518% 


162.736 


0.115% 


2.136 


0.261% 


22 


2.108464 


0.074% 


32.958 


0.526% 


164.967 


0.106% 


2.016 


0.095% 



an offset in long i tude a nd in latitude has actually been detected 
dTiscareno et al.l (120091) '). So, detection of non-hydrostatic con- 
tributions from observation of Mimas' orientation should not a 
priori be excluded. 



Appendix A: Notations used in the paper 
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7. Conclusion 

We have presented a theoretical study of the rotation of Mimas, 
in considering the 3 degrees of freedom of the rigid rotation, and 
different possible interior models, in assuming Mimas to be in 
hydrostatic equilibrium, or not. Moreover, we have considered a 
complete orbital motion, and also investigated the influence of 
tides on the rotation of Mimas. 

We estimate the physical longitudinal librations to have an 
amplitude of about 0.5°, i.e. nearly 2 km, the exact value depend- 
ing on the internal structure of Mimas. For a hydrostatic Mimas, 
a dense core lowers this amplitude. Non-hydrostatic contribu- 
tions are shown to be small as expected from Mimas shape in 
near hydrostatic equilibrium. Moreover, we expect an obliquity 
between 2 and 3 arcmin, while the polar motion can be ne- 
glected. The tidal deviation of Mimas' long axis should be neg- 
ligible as well, while this is the most inner main Saturnian satel- 
Ute. 

The Cassini spacecraft has already completed its initial four- 
year mission and the first extended mission, with a limited num- 
ber of Mimas flybys. Its orbit close to Saturn makes Mimas a 
diflicult target for Cassini observations. Since September 2010 
Cassini is in a second extended mission called the Cassini 
Solstice Mission during (and especially at the end) of which 
Cassini will likely have additional Mimas observations. We hope 
that future observations of Mimas will allow us to constrain its 
rotation and to get clues on its internal structure and orientation. 
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R, R, 
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m 


Mass of Mimas 


a> b > c 
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A<B<C 


Moments of inertia of Mimas 


A,. <Br<Cc 


Moments of inertia of the core of Mimas 


Pc > p., 


Densities of the core and the shell of Mimas 


kf 


Fluid Love number of Mimas 


ki 


Love number of Mimas 


Ji - -C20, C22 


Gravity coefficients of Mimas 


Proper modes 


A, u), (p, <I> 


Orbital modes 


n 


Orbital and spin frequency of Mimas 


U, V, w 


Rotational modes, respectively in longitude, latitude and wobble 


U)u, (JJv, OJw 


Frequencies of the rotational modes 


Rotation variables 


/, g, h, L, G, H 


Andoyer variables and the moments associated 


p, r, ^„ P, n, ij, 


Modified Andoyer variables and the moments associated 


Rotation outputs 




Tidal longitudinal librations 


r 


Physical longitudinal librations 


1 


Latitudinal librations 


e 


Obliquity of Mimas 


Qx + 1Q2 


Polar motion of Mimas 


Tidal data 


A 


Dissipative rate 




Displacement of the equilibrium 
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